Anomalous magnetoresistance on the topological surface 
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We investigate charge transport in two-dimensional ferromagnet /feromagnet junction on a topo- 
logical insulator. The conductance across the interface depends sensitively on the directions of the 
magnetizations of the two ferromagnets, showing anomalous behaviors compared with the conven- 
tional spin-valve. This stems from the way how the wavefunctions connect between both sides. 
It is found that the conductance depends strongly on the in-plane direction of the magnetization. 
Moreover, in sharp contrast to the conventional magnetoresistance effect, the conductance at the 
parallel configuration can be much smaller than that at the antiparallel configuration. 

PACS numbers: 73.43.Nq, 72.25.Dc, 85.75.-d 



Spintronics aims to manipulate and/or use the spin 
degrees of freedom in device functions. There are two 
mainstreams in spintronics: the control of charge trans- 
port by spins [H, H, B, B, [1, B , and the control of spins 
by the electric field @, [1,0, [lO, 11, HI- In the former, 
giant magnetoresistance [2| and tunnelling magnetoresis- 
tance [3, 4, 5] in metallic spin valves have received much 
attention [6|. In the latter, on the other hand, the spin- 
orbit interaction (SOI) plays an essential role to connect 
the charge and spin degrees of freedom. However, the 
role of the SOI in the magnetoresistance has not been 
considered seriously thus far. From this viewpoint, the 
recently discovered topological insulator offers an inter- 
esting laboratory to search for the possible spintronics 
functions with the strong SOI. 

Recent theoretical and experimental discovery of the 
two-dimensional quantum spin Hall system [3, 14. 15. 
E 0, 03, [Til . Ho], and its generalization to the topo- 
logical insulator in three dimensions [H, 21, 22|, 23. 24] 
have established the new state of matter in the time- 
reversal symmetric systems. The topological order in 
the bulk with the gap dictates that there should be 
the one-dimensional channels along the edge of the two- 
dimensional sample, or the two-dimensional metal on the 
surface of the three dimensional sample. These edge and 
surface states are protected by the time-reversal symme- 
try and the topology of the bulk gap, and are robust 
against the disorder scattering and electron-electron in- 
teractions. 

In topological surface state on 3D topological insula- 
tor, the electrons obey the 2D Dirac equations. This cor- 
responds to the infinite mass Rashba model [2j|, where 
only one of the spin-split bands exists. This has been 
beautifully demonstrated by the spin- and angle-resolved 
photoemission spectroscopy [2Q, [27]]. Therefore, the next 
step is to unveil the unique property of the surface state 
of the topological insulators, in particular that relevant 
to magnetism [28|, |2j| [lo|, HH, [32[ . One remarkable feature 
of the Dirac fermions is that the Zeeman field acts like 



vector potential: the Dirac Hamiltonian is transformed 
as k • a — > (k + H) • a by the Zeeman field H. There- 
fore, we can expect anomalous spin related property by 
the magnetic field in topological insulator. This clearly 
contrasts with the Dirac fermions on graphene since a is 
pseudospin there. Q 




FIG. 1: (Color online) schematics of F1/F2 junction. The fer- 
romagnetism is induced in the topological surface state due to 
the proximity effect by the ferromagnetic insulators deposited 
on the top. The current flows on the surface of the topological 
insulator. 

In this paper, we study charge transport in 2D topolog- 
ical ferromagnet /feromagnet junction. The ferromagnet 
is made of the topological surface with a ferromagnetic 
insulator on the top. We uncover anomalous magnetore- 
sistance in this spin- valve: the conductance strongly de- 
pends on the in-plane rotation with respect to the other 
magnetization direction. Moreover, in sharp contrast to 
the conventional magnetoresistance effect, the conduc- 
tance may have its minimum at the parallel configura- 
tion, while it may take a maximum near antiparallel con- 
figuration. This is due to the connectivity of the wave- 
function across the junction. 

We consider 2D ferromagnet /feromagnet junctions 
which is abbreviated as F1/F2 in the following. We focus 
on charge transport at the Fermi level inside the bulk gap 
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of the topological insulator, which is described by the 2D 
Dirac Hamiltonian 

H= f rn z k x +m x -i(k y + my)\,^ 

\ k x + m x + i(k y + m y ) -m z y 

where m x , m y and m z are exchange field and we set 
vp = h = 1. The ferromagnetism is induced due to 
the proximity effect by the ferromagnetic insulators de- 
posited on the top as shown in Fig. [H The interface is 
parallel to y-axis and located at x = 0. We choose the 
exchange field in the Fl side as mi = (m x ,m y ,m z ) = 
mi (sin cos sin 6 sin cp, cos 6) while in the F2 side, we 
set m x = rriy = and m z — m 2 . In actual experi- 
ment, one can use a magnet with very strong easy axis 



anisotropy for F2, and a soft magnet for Fl which can 
be controlled by a weak magnetic field. 

We consider the juction between different ferromag- 
nets. This type of interface should contain a built-in 
electric field. Thus, we take into account the potential 
drop V in F2 which represents the difference of the Fermi 
energies in the two ferromagnets. Also, due to the mis- 
match effect, some barrier region may be formed near 
the interface. We describe this region of the length L 
by the Dirac fermion with the barrier potentail U . Note 
that the potentials V and U may be tunable by gate elec- 
trode. Then, with the above Hamiltonian, wave function 
in the Fl side is given by 



ijj(x < 0) = 
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while the wave function in the barrier region is given by 

be 



i/j(0 < x < L) = ae ik " x I k " %ky 
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and that in the F2 side reads 



ip(x > L) 
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with E' 



V, where E 



y/m* + (k x + m x ) 2 + (k y + m y ) 2 = - ^fef + fe2 + 



U 



±^/m| + fei 2 + fc2 + V. Here, ± sign corresponds to the 

upper and lower bands. Below, "n" and "p" mean that 
the Fermi level crosses the upper and the lower bands, 
respectively. Also, r and t are reflection and transmission 
coefficients, respectively. It should be noted that the 
Fermi surface in the Fl is shifted by (—m x , —m y ) from 
the origin. Due to the translational invariance along 
the 7/-axis, the momentum k y is conserved. Hence, the 
common factor e lkyV is omitted above. 

By matching the wavefunctions at the interface x = 
and L, we obtain the transmission coefficient t. We con- 
sider the situation that the barrier region is sufficiently 
narrow so that we can take the limit of U — » oo and 
L — > while keeping Z = UL = const. Here, we omit 
the expression of t because it is rather complicated but 
we note that it contains the barrier parameter Z only in 
the form of cos Z and sin Z. Consequently, the transmis- 
sion probability and hence the conductance are 7r periodic 
with respect to Z. In the presence of Z, the spin direc- 
tion of wavefunction rotates through the barrier region, 
similar to the spin transister. [7] Thus, with increasing 



Z, the connectivity of the wavefunction changes, which 
crucially influences the conductance. 
^ We parametrize k x -\-m x = kp cos 0, k y -\-m y = kp sin <\>. 
Then, we have E = \Jm 2 z + (k x + m x ) 2 + (k y + m y ) 2 = 
y/m 2 z + k 2 F . 

Finally, we obtain the normalized tunneling conduc- 
tance as 



|t| 2 Re 



E 



(5) 



Now, let us discuss the applicability of our model. Typ- 
ical value of induced exchange field due to the magnetic 
proximity effect would be 5^50 meV 33l although 
this depends on the interface property and the material 
choice of the ferromagnet. On the other hand, E can be 
tuned by gate electrode or doping below the bulk energy 
gap (~ 100 meV). [35| Due to the presence of the ferro- 
magnet, time reversal symmetry is broken. This would 
tame the robustness against disoder. However, high qual- 
ity topological insulator can be fabricated now, the mean 
free path of which is sufficiently large, |3(| and hence lo- 
calization does not occur in the temperature region of 
our interest and surface state is stable for exchange field 
smaller than the bulk energy gap. The above approxima- 
tion of discontinuous change of potential can be justified 
as follows. The characteristic length of the wavefuction 
is £ = hvp I m z while Thomas- Fermi screening length A is 
given by 1/A = e 2 N(E) in 2D where N(E) is the density 
of states at the Fermi level. Then, we have £/A ~ E/m z 
using vp — 6 x 10 5 m/<s for Bi2Se3 [35]. Thus, we obtain 
£ > A for E > m z . 

In the following, we will show results for mi = \^X9E. 
The tunneling conductance strongly depends on how the 
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FIG. 2: (Color) tunneling conductance a with Z = for 
TU2 — ((a) and (c)), and mi — V0.9E ((b) and(d)). n-n 
junction at V = in (a) and (b). p-n junction at V = 2E in 
(c) and (d). 
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FIG. 3: (Color) similar plots to Fig. [2] with Z = tt/2 for 
m2 = ((a) and (c)), and 1712 — V&9E ((b) and(d)). n-n 
junction at V — —E in (a) and (b). p-n junction at V — 2E 
in (c) and (d). 

wavefunctions connect between both side, which we will 
explain with Fig. [2] for Z = and Fig. [5] for Z = tt/2. 
To describe the physics, we first consider the Z = case. 

In Fig. [2j we show the normalized tunneling con- 
ductance (7 in n-n junction for (a) Vfi2 = and (b) 
m2 = V0S)E. In Fig. [2] (a), the F2 is no more ferromag- 
netic. Nevertheless, the conductance strongly depends on 
the direction of the magnetization in the Fl. At = or 
7r, the mismatch of the wavefunctions between the two 
sides and that of the sizes of Fermi surfaces suppresses 
a, because the energy E is near the bottom of the upper 
band in Fl while there is no gap in F2. At 6 = 
on the other hand, the wavefunctions and the sizes of 
the Fermi surfaces are the same on both sides except 
the shift of Fermi surface in the momentum space due to 



the in-plane component of the magnetization as shown in 
Fig. [4j However, this misfit of the in-plane momentum 
between the two sides gives rise to a strong dependence 
of (7 on the in-plane rotation angle (/?, which is not seen 
in the conventional magnetoresistance effect. Since k y 
is conserved, the positions of the Fermi surfaces strongly 
influence the charge transport: if exchange field points to 
x-axis, there is no evanescent wave. On the other hand, 
when exchange field is applied in ^-direction, the Fermi 
surface moves to the k y direction and hence the overlap 
region between fc^'s in the Fl and F2 is reduced. There- 
fore, the number of the evanescent modes increases and 
hence the conductance is strongly suppressed. Thus, we 
can obtain giant magnetoresistance in this system. 

In Fig. [2] (b), the conductance is large at the paral- 
lel configuration (6 = 0) while it is small for antiparallel 
configuration [6 = tt) . This dependence, similar to the 
conventional magnetoresistance effect 0, Q, can be un- 
derstood by the overlap intergral of the wavefunctions on 
both sides, as discussed later. Note that at the antipar- 
allel configuration, the domain wall structure generates 
the edge state at the interface. However, this edge state 
is merged with the surface state and do not make a sig- 
nificant contribution to the conductance. 

We show tunneling conductance in p-n junction with 

V = 2E for m 2 = in FigE] (c) and m 2 in 
Fig. [2] (d). In Fig. [2] (c), a similar tendency to Figure 
Efa) is seen. In Fig. [2](d), in stark contrast to the con- 
ventional magnetoresistance effect, the conductance takes 
minimum at the parallel configuration (6 = 0) while it 
takes maximum near antiparallel configuration (6 = tt). 

To understand these results intuitively, we describe the 
underlying physics in Fig. [5] where the arrows indicate 
the spin directions in the limiting case of \m z \ — > oo, 
showing the connection of the wavefunctions on both 
sides. In the n-n junctions, the tunneling amplitude is 
determined by the overlap of the same eigenfunctions for 
parallel configuration (6 = 0), while for antiparallel con- 
figuration (6 = tt) it is given by the overlap of the differ- 
ent eigenfunctions, as shown in Figs. [5] (a) and (b). Thus, 
the tunneling amplitude takes its maximum at = 0, and 
this explains the dependence of the conductance in Fig. 
[2] (b). In a similar way, in p-n junctions, we find that the 
tunneling amplitude becomes larger at = tt than that 
at = as shown in Figs. [5] (c) and (d). This is the ori- 
gin of the anomalous dependence of the conductance in 
Fig. El(d). 

Next, we consider the influence of the barrier poten- 
tial Z and the potential drop V. Figure [3] exhibits tun- 
neling conductance with Z = tt/2 for n-n junction at 

V = — E ((a) and (b)), and p-n junction at V = 2E 
((c) and (d)), which should be compared with Figs. [2] 
(a), (b) and (c), (d), respectively. In the n-n junction, 
at V = — E the Fermi surface becomes larger than that 
at V = 0. Then, the effect of the shift of the Fermi 
surface becomes less important. This leads to the weak 
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FIG. 4: (Color online) Positions of Fermi surfaces. On the Fl, 
the Fermi surface moves as illustrated, as mi rotates around 
z-axis. 




FIG. 5: (Color online) Connectivity of the wavefunction 
across the n-n junction ((a) and (b)), and p-n junction ((c) 
and (d)). The magnetizations are parallel (0 = 0) in (a) and 
(c), while they are antiparallel (0 = ty) in (b) and (d). The 
arrows represent the electron's spin. 

if dependence as shown in Figs. [3] (a) and (b). Note 
that the spin of the eigenstate of Eq.(l) is parallel to 
(k x + m x , k y + m y , m z y. Therefore, the in-plane compo- 
nent of the spin in the wavefunction in F2 is dominant for 
V = — E. The conductance is largest at = 7r/2, when 
the in-plane spin component in Fl and hence the overlap 
of the wavefunctions between Fl and F2 are maximum 
(Figs. [3] (a) and (b)). In Fig. [3] (c), since ra 2 = 0, the 
spin rotation by Z does not make a significant change in 
the conductance compared to Fig. [2](c). At Z = tt/2 the 
spin is half rotated, and therefore the tendency becomes 
opposite comparing Fig. [2] (d) and Fig. [3] (d). 

In summary, we studied charge transport in 2D topo- 
logical ferromagnet /feromagnet junction. The ferromag- 
net is made of the topological surface with a ferromag- 
netic insulator on the top. We found anomalous magne- 
toresistance in this topological spin-valve. 
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